Long-time electron spin storage via dynamical suppression of 
hyperfine-induced decoherence in a quantum dot 
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The coherence time of an electron spin decohered by the nuclear spin environment in a quantum 
dot can be substantially increased by subjecting the electron to suitable dynamical decoupling 
sequences. We analyze the performance of high-level decoupling protocols by using a combination of 
analytical and exact numerical methods, and by paying special attention to the regimes of large inter- 
pulse delays and long-time dynamics, which are outside the reach of standard average Hamiltonian 
theory descriptions. We demonstrate that dynamical decoupling can remain efficient far beyond its 
formal domain of applicability, and find that a protocol exploiting concatenated design provides best 
performance for this system in the relevant parameter range. In situations where the initial electron 
state is known, protocols able to completely freeze decoherence at long times are constructed and 
characterized. The impact of system and control non-idealities is also assessed, including the effect 
of intra-bath dipolar interaction, magnetic field bias and bath polarization, as well as systematic 
pulse imperfections. While small bias field and small bath polarization degrade the decoupling 
fidelity, enhanced performance and temporal modulation result from strong applied fields and high 
polarizations. Overall, we find that if the relative errors of the control parameters do not exceed 
5%, decoupling protocols can still prolong the coherence time by up to two orders of magnitude. 



PACS numbers: 03.67.Pp, 76.60.Lz, 03.65.Yz, 03.67.Lx 



I. INTRODUCTION 

Electron and nuclear spin degrees of freedom are 
promising candidates for a variety of quantum informa- 
tion processing (QIP) devicesii^. While the wide range 
of existing microfabrication techniques make solid-state 
architectures extremely appealing in terms of large-scale 
integration, such an advantage is seriously hampered by 
the noisy environments which are typical of solid-state 
systems, and are responsible for unwanted rapid deco- 
herence. For an electron spin localized in a GaAs quan- 
tum dot (QD)^, for instance, the relevant coherence time 
is extremely short: at typical experimental temperatures 
(T ~ 100 mK) and sub-Tesla magnetic fields, the elec- 
tron free induction decay (FID) time T2 ~ 10 ns^i^ii, 
the dominant decoherence mechanism being the hyper- 
fine coupling to the surrounding bath of Ga and As nu- 
clear spins. Although efforts for achieving faster gating 
times may contribute to alleviate the problem, it remains 
both highly desirable and presently more practical to ex- 
tend the coherence time of the central system (the elec- 
tron spin) in the presence of the spin bath. 

Several proposals have been recently put forward to 
meet this challenge. A first strategy is to manipulate the 
spin bath. Polarizing the nuclear spins, for instance, may 
significantly increase the coherence time^iS, provided that 
nuclear-spin polarization > 99% may be achieved. This, 
however, remains well beyond the current experimen- 
tal capabilities. Another suggestive possibility, based on 
narrowing the distribution of nuclear spin state o^"'^^'^^ , 
has been predicted to enhance electron coherence by up 
to a factor of hundred, upon repeatedly measuring and 



pumping the electron into an auxiliary excited trion state 
- which also appears very challenging at present. 

As an alternative approach, direct manipulation of the 
central spin by means of electron spin resonance (ESR)^^ 
and dynamical decoupling (DD)^ii^ii^iii techniques ap- 
pears ideally suited to hyperfine-induced decoherence 
suppression, in view of the long correlation time and non- 
Markovian behavior which distinguish the nuclear spin 
reservoir. A single-pulse Hahn-echo protocol has been 
implemented recently in a double-QD device^, increasing 
the coherence time by two orders of magnitude. Signifi- 
cant potential of more elaborated pulse sequences, such 
as the multi-pulse Carr-Purcell-Meiboom-Gill (CPMG) 
protocol's and concatenated DD {CVtDp.-^^^^^ has 
been established theoretically for a single QD subjected 
to a strong external bias field, whereby the electron ef- 
fectively undergoes a purely dephasing process. The DD 
problem for the more complex situation of a zero or low 
bias fields, where pure dephasing and relaxation compete, 
has been recently examined in Refs. Having es- 

tablished the existence of highly effective DD schemes for 
electron spin storage, the purpose of this work is twofold: 
first, to gain a deeper understanding of the factors infiu- 
encing DD performance and the range of applicability 
of conclusions based on analytical average Hamiltonian 
theory (AHT) approaches; and second, to assess the in- 
fluence of various factors which may cause the system 
and/or control Hamiltonians to differ from the idealized 
starting point chosen for analysis. 

Aside from its prospective practical significance, devel- 
oping and benchmarking strategies for decoherence sup- 
pression in various spin nanosystems is interesting from 
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the broader perspective of quantum control theory. In 
particular, standard theoretical tools usually employed 
for the analysis of DD performance, such as AHT and the 
Magnus Expansion (ME), have very restrictive formal re- 
quirements of applicability (very fast control time scales, 
bounded environments, etc), which may be hard to meet 
in realistic systems. Thus, in-depth studies of physically 
motivated examples are essential to understand how to 
go beyond the formal error bounds sufficient for conver- 
gence, and to identify more realistic necessary criteria 
for DD efficiency. In this sense, a QD system, described 
by the central spin model (a central spin-1/2 interact- 
ing with a bath of N external spin a^^'^^ ) both provides a 
natural testbed for detailed DD analysis, and paves the 
way to understanding more complex many-spin central 
systems. 

In this work, we present a quantitative investigation 
of DD as a strategy for robust long-time electron spin 
storage in a QD. The content of the paper is organized 
as follows. In Sec. |TT1 we lay out the relevant control 
setting, by describing the underlying QD model as well 
as the deterministic and randomized DD protocols under 
consideration. Among periodic schemes, special empha- 
sis is devoted to the concatenated protocol (PCDD2), 
which was identified as the best performer for this sys- 
tem in Ref. [H. Exact AHT results are obtained up to 
second-order corrections in the cycle time, which allow 
additional physical insight on the underlying averaging 
and on DD-induced bath renormalization to be gained. 
Details on the methodologies followed to assess the qual- 
ity of DD and to effect exact numerical simulations of the 
central spin coupled to up to Af = 25 bath spins are also 
included in Sec. |TT1 

In Sec. mil numerical results on best- and worst-case 
performance of DD protocols for short evolution times are 
presented and compared with analytical predictions from 
AHT/ME under convergence conditions - in particular, 
LOcT <C 1, where ujc and r are the upper cutoff frequency 
of the total system-plus-bath spectrum, and the time in- 
terval between (nearly instantaneous) consecutive con- 
trol operations, respectively. Evolution times as long as 
~ IOOOT2 for T < T2* are able to be investigated numer- 
ically, other decoherence mechanisms becoming relevant 
for yet longer times. In the best-case scenario, where de- 
coherence of a known initial state may be frozen under 
appropriate cyclic DD protocola^^, the dependence of the 
attainable asymptotic coherence value on r is elucidated 
in the small r region. Sec. IIVI is devoted to further in- 
vestigating the effect of experimentally relevant system 
features and/or non-idealities, such as the presence of 
residual dipolar couplings between the nuclear spins, the 
influence of an applied bias magnetic field, and the role of 
initial bath polarization. In Sec. |Vl the effect of system- 
atic control imperfections such as finite width of pulses 
and rotation angle errors is quantitatively assessed. We 
present conclusions in Sec. I VII 



II. SYSTEM AND CONTROL ASSUMPTIONS 

In this section we describe the model spin-Hamiltonian 
for a typical semiconductor QD, and typical picture of the 
decoherence dynamics of an electron spin in a QD. We 
present the DD methods used to suppress the electron 
spin decoherence and the metrics for DD performance, 
followed by numerical methods we employed. 

A. QD model Hamiltonian 

The decoherence dynamics of an electron spin S local- 
ized in a QD and coupled to a mesoscopic bath consisting 
of A^ nuclear spins Ifc, fc = 1, . . . A^, may be accurately de- 
scribed by the effective spin Hamiltonian— 

H ^Hs + Hb+Hsb, (1) 

where, in units h = 1, 

Hs = u;aS,, (2) 

N N 

Hb =Y.T. rfc£(Ifc • h - 31^1!), (3) 
k=i e<k=i 

N 

HsB = ^AfeS-Ife, (4) 

k=l 

Here, the Hamiltonian Hs describes the electron spin, 
LjQ = g*fiBBQ being the Zeeman splitting in an exter- 
nal magnetic field Bq, g* is the effective Lande fac- 
tor of the electron spin, and /xs is the Bohr's magne- 
ton. The Hamiltonian Hb is the bath Hamiltonian, de- 
scribing dipolar interactions with strength Tki between 
nuclear spin k and £. For our analysis, a nearest- 
neighbor approximation to Hb is adequate, with Tki 
taken as random numbers uniformly distributed between 
[— ro,ro], so that the parameter Fq upper-bounds the 
strength of intra-bath couplings. The relevant system- 
bath coupling Hamiltonian Hsb accounts for the Fermi 
contact hyperfine interaction, the coupling parameter 
Ak = (87r/3)(7|/iB5„/i„M^(xfc) being determined by the 
electron density u^(xfc) at the fc-th nuclear spin site 
and by the Lande factor of the nuclei Other small 
contributions to the total QD Hamiltonian, such as the 
Zeeman splitting of the nuclear spins and anisotropic 
electron-nuclear couplings, may be neglected for the cur- 
rent purposes^Si^ (see, however, Ref. |3l| for recent devel- 
opments on controllability in the presence of anisotropic 
hyperfine couplings). 

Upon tracing over the nuclear-spin reservoir, the elec- 
tron spin described by Eq. ([T]) undergoes fast decoherence 
with a characteristic FID time of 

" ~ A4/(/+l)' V A^ ' 

where A sa lO^^'/ieV for typical GaAs QDs with A^ = 
10^ Ij. = I = 3/2 nuclear spins^^'^^. This results in a 
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Tj* value of about 10 ns^i^ii^, which is too short for 
QIP apphcations. In simulations, we shall neglect the 
/-dependence of the FID time and simply set the nuclear 
spin value to Ik = 1/2 for all k. The FID time also 
turns out to depend very weakly on the applied bias field 
Bq as long as Bq is smaller than or comparable to the 
Overhauser field of the unpolarized nuclear spins^^. A 
similar conclusion holds for a weakly polarized nuclear 
spin bath, with polarization p <C 1. Throughout this 
paper, energies shall be expressed in units of A, and time 
shall be expressed in units of 1/A. Since Fq ^ A^. for 
typical GaAs QDs^, we shall take Fg = (hence Hb = 
0) unless explicitly stated. 



B. QD decoupling protocols 

Compared to other decoherence control strategies, DD 
has many attractive features: it is a purely open-loop 
control method which, as such, avoids the need of mea- 
surements and/or feedback; it does not rely on a par- 
ticular initial state which might be hard to prepare; its 
design and implementation may significantly benefit from 
the extensive expertise available from nuclear magnetic 
resonance (NMR) and ESR techniques. 

We shall first assume to have access to ideal control re- 
sources, and defer the discussion of control limitations to 
Sec. |Vl In this idealized setting, DD of the electron spin 
is implemented by subjecting it to a sequence of bang- 
bang TTfi pulses^^, each instantaneously rotating the spin 
S along a given control axis h by an angle tt. Consecu- 
tive pulses are separated by a time interval r, resulting 
in a total time-dependent Hamiltonian of the form 



Htot{t) = H,{t)+H, 



(5) 



where 5{-) is the delta function, £ ~ 1,2, .. . labels the 
applied pulses, and H is given by Eq. (jT]). The evolution 
of the coupled electron-nuclear system in the physical 
frame is then described by the unitary propagator 



i I Htot{s)ds 



U{t) = r exp 



where T denotes as usual temporal ordering. 



(6) 



1. Deterministic single-axis DD 

DD protocols involving control rotations of the central 
spin about a fixed axis achieve selective averaging of spin 
components perpendicular to the rotation axis. Such se- 
quences are effective when a preferred direction is present 
either in the underlying physical Hamiltonian or in the 
initial quantum state to preserve. For instance, in the 
presence of a strong static bias field, i?o ^ 1 T, electron 



spin flips are energetically suppressed, and pure dephas- 
ing in the transverse direction is the dominant decoher- 
ence sourcei2iS2ii2^iM. In such a situation, the interaction 
Hamiltonian Eq. ^ simplifies to 



Hcff — S; 



k 



(7) 



with renormalized coupling constants A'f.. The z- 
component of the central spin remains constant, while 
the transverse component s± — ^ (Sx)^ + (Sy)^ under- 
goes Gaussian FID with time constant By applying 
a TTj; pulse at times t = t and 2r (the so-called Hahn 
echo, for brevity denoted as [rXrX]), the electron spin 
is refocused completeljJ^ at time T = 2t, 



(8) 



A protocol consisting of repeated spin-echoes shall be 
referred to as CPMG henceforth. In the case of zero 
or small bias field where the full Fermi contact interac- 
tion is relevant (see Sec. IIV Cl for more discussion on the 
case of a nonzero bias field), selective DD no longer re- 
moves the effect of the hyperflne coupling on a generic 
initial electron state. Rather, to lowest order in t, the 
control has the effect of symmetrizing the system-bath 
Hamiltonian according to the direction of the applied 
pulsesiii^. The presence of such a control-induced ap- 
proximate symmetrization is essential to understand the 
possibility of decoherence freezing2^, to be further dis- 
cussed in Sec. IIII Al as well as in the Appendix. 



2. Deterministic two-axis DD 

In the absence of a bias field {Bq = 0), relaxation 
in the longitudinal direction is as important as dephas- 
ing. Thus, removing the effect of the nuclear reservoir 
on the electron spin is only possible by using a con- 
trol protocol which achieves non-selective (or universal) 
DD. Several sequences have been constructed for finite- 
dimensional systems by assuming control over a basic set 
of unitary operations forming a discrete group G = {gj}, 
j = 0, 1, . . . , 1^1 - 1 (DD group) i^iiiiia. In the simplest 
case of cyclic DD, the control propagator is sequentially 
steered through the DD group in a predetermined or- 
der, the change from gi to gj being effected through the 
application of a bang-bang pulse Pij — gjgj - for in- 
stance, G = {I, X} in the above-mentioned CPMG pro- 
tocol, / denoting the identity operation. Thanks to the 
inherent periodicity of the control action, with cycle time 
Tc = \G\t, the DD analysis can be naturally carried out 
within the AHT by invoking the MBHiiiO^, 



U{Tc) = exp{-inT^ 



n 



fc=0 



(k) 



(9) 
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where H. denotes the AHT. In principle, arbitrary high- 
order contributions TL^ ^ may be expUcitly evaluated by 
knowing the applied control sequence'^^. A DD scheme 
for an open-system Hamiltonian of the form ([1]) is said 

to be of order m if Ti!"^^ = Hb and the first non-zero 
contribution to Ti. arises from Ti'™^ , thereby being of or- 
der ©(r™) in the expansion parameter Tc- Estimates of 
the convergence radius for the ME depend sensitively on 
how the strength of H is quantified, which may be es- 
pecially delicate for mesoscopic- and infinite-dimensional 
environments^flii^i^. A conservative convergence bound 
arises by assuming that a finite upper spectral cutoff 
may be identified, \\H\\ ujc < oo, and the condition 
uJcTc < 1 is obeyed^Si^i^iiS. A more precise character- 
ization of sufficient convergence conditions for the ME 
has been recently established in Ref. 0. 

Periodic DD (PDD) is the simplest non-selective cyclic 
protocol, which ensures that the unwanted evolution is 
removed to first order in the ME at every T„ = nTc, n € 
N, for sufficiently short Tc. For a single central spin, PDD 
is based on the irreducible Pauli group Gp — {I, AT, Y, Z}, 
where X, Y, Z are Pauli matrices, up to irrelevant phase 
factorsi^iii. A possible implementation, corresponding 
to the group path {I,X,Z,Y), involves two-axis control 
sequences of the form 



Ci = PDD 



Cq X CqZCq X Co z 



Ooe Oo 
(10) 



with Co = exp(— ii?r) representing the operator of free 
evolution between pulses. Note that in the second line 
of Eq. (fTO]) , the time convention used for pulse sequence 
descriptions is followed, that is, time evolves from left to 
right. For such a sequence, the lowest-order ME reads 



exp 



-iTcin 



(0) 



(11) 



with 



n^°^ = ^[H^ + H3 + H2 + Hi]t = Hb, (12) 

c 



3=1 i=l 



= -^vlH E ^Miki! + iki!), (13) 

k e 

where Hi ^ H, H2 = XHX, H3 = YHY, and = 
ZHZ, respectively. Thus, PDD achieves first-order DD, 
with the unwanted hyperfine coupling vanishing in the 

limit T — > 0, and Ti!'^'^ = for a non-dynamical bath. 
Note that the decoupling happens after completion of 
each cycle, so everywhere in the subsequent analysis of 
cyclic DD, we shall assume that the evolution is sampled 
stroboscopically at instants tj = jTc,j = 1, 2, . . . 

A simple strategy to improve over PDD averaging is 
based on exploiting time-symmetrizationi^iii - leading 



to so-called symmetric DD (SDD). Corresponding to the 
above-defined PDD protocol, the SDD protocol relevant 
to our problem is defined by the control cycle 



SDD — CqXCqZCqXCqCqXCqZCqXCq . 



(14) 



Symmetrization guarantees that the electron spin opera- 
tors are cancelled in all the odd terms of the ME {k odd), 
at the expense of making Tc twice as long as in PDD. As 
long as the bath is non-dynamical, Hb = 0, this also 

yields Ti.^^^ — 0. As one may verify by direct calculation, 
the lowest-order term containing the electron spin is 



7(2) 



^3 8 fe J 



6T E E E( [^^-' t^^' ^*]] + 
^ k=\ j=l 1=1 



(15) 



A powerful method to enhance DD performance by 
systematically shrinking the norm of higher-order terms 
is to resort to concatenated designi^. CDD relies on a 
temporal recursive structure, so that at level i+\ the pro- 
tocol is defined recursively as Cf+i = CiXCiZCiXC^Z, 
Cq denoting as before free evolution under H . Building 
on the results of Ref. '2^, we focus in this work on a trun- 
cated version of CDD, whereby concatenation is stopped 
at a certain level and the resulting supercycle is repeated 
periodically afterwards. Such a protocol is referred to as 
PCDDf, with Tc = 4V (note that i=l recovers PDD). 

Specifically, a cycle of the PCDD2 protocol consists of 
two identical half-cycles, and has a form 

[CoXCoZCoXCoYCoXCoZCoXCo][repeatl (16) 

where [•] denotes each half-cycle. This protocol leads 
to especially remarkable DD performance and averaging 
structure for the central-spin problem. First, because the 
protocol is time-symmetric, averaging is at least of second 

order, as in SDD that is, n^°^ = 0, and h''^^ = 0. How- 
ever, in contrast to SDD, the second-order contribution 

(2) 

H turns out to be an effective pure bath term which 
renormalizes the bath Hamiltonian without directly con- 
tributing to the decoherence dynamics. Thus, operators 
mixing electron and nuclear spins can only appear at or- 

— (4) 

der Ti. and higher. The fact that a pure-bath contri- 
— (2) 

bution H is generated to second-order in the ME, and 
that a particularly favorable scenario is to be expected 
for CDD convergence, was implied in Ref. I20!. Remark- 
ably, two additional features emerge to second order in 
the controlled dynamics: 

(i) The DD-induced bath Hamiltonian has a regular 
coupling structure. 



2 16 fe j 

--^ E E E ( t^'^' t^^-' ^']] + [[^fc' ^^i' 

k=l j=l 1=1 
N N 



=E E r'.,(i.-i, -3/f/;), 



(17) 



2—1 j<ci—l 
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where 

i ,j = T ^ . (18) 

That is, to second order in r, the hyperfine interaction 
between the electron and the nuclei is removed, and an 
effective dipolar Hamiltonian with control-renormalized 
couplings is induced on the nuclear spins ~ compare with 
the Hamiltonian in Eq. 

(ii) To the same level of accuracy in r, it is possible 
to show that a similar averaging is achieved by a half- 
PCCD2 protocol - that is, a protocol whose cycle consists 
of just the first half in ([T6|) . In practice, this may be 
useful to reduce the number of required pulses for a given 
desired accuracy. 

3. Randomized DD 

Randomized design offers another approach to improve 
DD performance, by both ensuring robust behavior in the 
presence of intrinsically time- varying open-system Hamil- 
tonians and by minimizing the impact of coherent error 
accumulation at long evolution times^i^i^i^. Unlike 
deterministic schemes, randomized DD is distinguished 
by the fact that the future control path is not fully known 
in advance. Analysis is most conveniently carried out in 
a logical frame that follows the applied control^", and by 
tracking the applied sequence to ensure that appropriate 
frame compensation can be used to infer the evolution in 
the physical frame. Although loss of periodicity in the 
control Hamiltonian causes AHT not to be directly appli- 
cable, contributions of different orders in r may still be 
identified in the effective Hamiltonian describing a given 
control sequence. 

Among all randomized protocols^^iiid^, we consider 
a few representative choices. Naive random DD (NRD) 
corresponds to changing the control propagator accord- 
ing to a path which is uniformly random with respect to 
the invariant Haar measure on Q. For our system, this 
means that a, X, Y , Z pulse, or a no-pulse is applied with 
equal probability at every instant tj — jr, j = 0, 1, . . . . 
So-called random path DD (RPD) merges features from 
both pure random and cyclic design, by involving FDD 
cycles where the path to traverse Q is each time chosen 
at random among the \Q\\ possible ones. For our sys- 
tem, we implement a simplified pseudo-RFD protocol by 
restricting to cycles which always begin with the iden- 
tity - that is, at every instant r„ = 4nr, n = 1,2,..., 
we randomly choose a sequence from the following set of 
{\Q\ — 1)! possibilities: 

CqXCqYCqXCqY, CoYCqXCqYCqX, 
CoYCqZCqYCqZ, CqZCqYCoZCqY}. 

Unlike NRD, RFD guarantees that averaging of Hsb 
is retained to lowest order over each control cycle. By 



enforcing symmetrization on the RPD protocol, sym- 
metrized random path DD (SRFD) is obtained, whereby 
every sequence from the above set is augmented by its 
time-reversed counterpart. SRFD additionally ensures 
that all odd terms in the effective Hamiltonian which in- 
volve the electron spin operators disappear at Tp = Spr, 
p=l,2,.... 



C. DD performance metrics 

In order to meaningfully compare different DD proto- 
cols for given control resources (in particular, finite t), an 
appropriate control metric should be identified. Different 
choices may best suit different control scenarios. Quan- 
tities such as pure-state input-output fidelity and purity, 
for instance, depend strongly on the initial state of the 
central system, thus being appropriate when preservation 
of a known electron initial state is the intended DD goal. 
Average input-output fidelity and gate entanglement fi- 
delity do not rely on a particular initial state, however 
they depend on the probability distribution of the initial 
states and quantify typical performanc o'^^'^^i'^^ . Taking 
advantage of the small dimensionality of the central spin 
system, an accurate way for quantifying DD performance 
at preserving an arbitrary electron pure state is to evalu- 
ate both best-case and worst-case input-output fidelities. 

For a given initial state \4>s{0)), recall that the input- 
output fidelity is defined as 

F{t) = Tr[ps{t)ps{0)] - (V'sWIpsWI^sW), (19) 

where ps(t) is the reduced density matrix of the central 
spin S at time t, ps{t) — TTBp{t), p(t) and Trs denoting 
the total density operator at time t, and the partial trace 
over the bath degrees of freedom, respectively. F{t) gives 
a measure of how far the central system has evolved away 
from its initial state. The best-case (6) and the worst- 
case (w) fidelities are then naturally defined as 

Fb{t) = max|^^){F(i)}, F^{t) = min|^^){F(i)}. 

D. Numerical methodology 

Analytical bounds on the expected worst-case fi- 
delity decay for various DD protocols have been ob- 
tained for sufficiently short evolution times based on 
either AHT and/or on additional simplifying assump- 
tion o"'^^'^°i'*°'^'^i'^^'^^ . For cyclic schemes, the relevant error 
bounds assume convergence of the ME series, thus requir- 
ing, in particular, that tUcTc < 1, where in our problem 
ujc is the highest frequency of the total (electron plus nu- 
clei) spectrum. For a typical GaAs QD with iV ^ 10^ 
bath spins, 

...V^.i^^20GHz, 

k 
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thus strict convergence of the ME imphes extremely short 
characteristic control time scales, r ~ 10 ps. Moreover, 
for ESR experiments under resonance conditions, ojc ^ 
20 GHz is about 400 times larger than currently available 
carrier pulse frequency^^, and roughly 20 times larger 
than attainable exchange-gating frequency^^. 

In order to evade the strict convergence requirement 
and to study DD performance in regimes of more direct 
experimental relevance, numerical simulations are neces- 
sary. Specifically, we are interested in pulse separations 
of the order of l/cr, where 




is the characteristic half-width of the coupling spectrum, 
as opposed to the highest frequency cjc, with a being 
roughly \/N ^ 10'^ times smaller than ujc- Thus, > 4t 
~ ^/Nuj~'^ , with typical inter-pulse delays values of the 
order of 10 ns, which is not too far (within an order 
of magnitude) from current experimental capabilities*^. 
Furthermore, in order to access DD in the long-time limit, 
we consider up to thousands of control cycles. Nuclear 
spin environments consisting of iV < 25 spin- 1/2 are in- 
vestigated (with a corresponding Hilbert space dimension 
of ~ 7 - 10'') - giving hope that our main conclusions may 
be extrapolated to real mesoscopic environments. 

The initial state of the entire system is taken to 
be a direct product of the electron initial spin state 
Ps{0) = |'0s(O))('0s(O)l and the bath initial spin state 
Pb(0). In most cases, we assume an unpolarized spin 
bath, described by the thermal equilibrium density op- 
erator pb(0) = 2~^/2ivx2«i where /2™x2™ is the 2^- 
dimensional identity matrix (polarized initial bath states 
shall be considered in Sec. IIV Cp . Such a maximally 
mixed spin state reflects the fact that for typical experi- 
mental dilution-refrigerator temperatures (T ~ 100 mK), 
the thermal energy is much larger than the energy scale 
of the intra-bath spin interactions. For a small number of 
bath spins {N < 8), we directly simulate the evolution of 
the total system's density operator, followed by a partial 
trace over B. For larger N, this approach is not computa- 
tionally efficient, and we perform simulations by assum- 
ing that the total system is in a pure stat o^'^i^'^'^^i^^ . In 
this case, psifi) is approximated as a uniformly random 
superposition of environment product states, |V'b(0)) = 

Si=iCi|'/'i) of all possible tensor products of the form 
|<^) = I T)i (X" I 1)2 <8i I Da <Xi ■ • ■ <Xi I T)Ar, where a are uni- 
formly distributed random number a^'^i^^'^^i^^'^^ , subject 
to |cip = 1. Thanks to concentration of measure on 
the space of random pure states, such an approximation 
has an exponentially small error, of the order of 2~^/^, 
with respect to using the identity state - that is, ~ 0.5% 
for N = 15. 

The desired best /worst case fidelity F^,/^ is evaluated 
by invoking quantum process tomography^'^. Four differ- 



ent initial states of the electron are considered, 

m = I T), 

m = 

1^.) = ^(|T) + U)), 
1^.) = ^(IT)+*U)) 

and for each such state the time-dependent Schrodinger 
equation of the total system subjected to DD is solved, by 
employing the Chebyshev polynomial expansion method 
to calculate the evolution operato r^^'^^ . At the final evo- 
lution time T, the four reduced density matrices (ob- 
tained upon partial trace over B) are used to compute 
the superoperator matrix Xmn{T) which describes the 
electron spin dynamics according to the linear map 

3 3 

m— n— 

where Kq = /, i^i = X, K2 = -iY, Ks = Z, and 
Ps specifies evolution in the logical frame, with j55(0) = 
P5(0). Since an arbitrary initial pure state of S on the 
Bloch sphere may be parameterized as 

1^-5(0)) -cos(0/2) I T)+sin(0/2)e'^| i), 

with 9 e [0,7r] and ip € [0, 27r], Fij/^{T) simplifies to 

F^,^{T) = max/min(,^^) (V's(0)|ps(T)|^5(0)) . 

In practice, after determining the matrix Xi we obtain 
Fb/w{T) by numerical optimization, using a statistically 
meaningful set of initial guesses for (0, 1^9). 



III. RESULTS ON DD PERFORMANCE 

In this section we summarize results on best- and 
worst-case performance of deterministic DD protocols 
relevant to the QD problem, including analytical esti- 
mates at short times. The discussion of randomized pro- 
tocols is postponed to Sec. IIIIBl as there is little differ- 
ence between best and worst case. 



A. Best-case performance and decoherence freezing 

Investigation of best-case performance both provides 
an upper bound on the attainable DD fidelity, and reveals 
a remarkable phenomenon of single-qubit deterministic 
DD: the optimal performance bound is reached for initial 
states of the electron, which depend on the geometry of 
the control protocol and cause the long-time fidelity to 
freeze at a non-zero value determined by the DD rate^^*^. 
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1. Short-time behavior 

In order to bridge analytical and numerical results, we 
begin by considering the performance of a single DD cy- 
cle, in the limit where the inter-pulse delay r is small 
enough to justify the application of AHT and the ME. 
Clearly, perfect (stroboscopic) preservation would be en- 
sured for an initial electron spin state which is an eigen- 
state of the full AHT Ti. For finite DD accuracy, the best- 
case scenario still corresponds to initializing the electron 
in an eigenstate of the lowest-order term of the ME, so 
that the fidelity decay after one cycle is determined by 
the next-lowest-order contribution to Ti.. 

Let the total cycle propagator be expanded as 



UiTc) = exp[-zT,(H^*'^ + H*'^ + W*'^ 



•)], (21) 



and consider, to begin with, the simplest [CqZCqZ] 

CPMG protocol. Then H^"^ ^ (g> Y^k^^Ik, and the 

next dominant contribution arises from Ti}^\ By using 
Eq. (jl9p . and by assuming that |V's(0)) is invariant under 



Fb{T,) 



for best-case performance, one finds 

CPMG _ 1 rr.2 



Pb{0) [aH^'^ 



1 - T^Tyb 
1- aT^T^ + 0{T^T^), 



\i's{0)) 



(22) 



where 



\^s) 



(^5|(H^'^)2|^S) 



denotes the partial variance of the dominant correction 

Ti in the initial electron state iV's), and ar^, a G M, 
gives the expectation of such variance operator over the 
initial bath state. Thus, a is clearly protocol- and bath- 
state dependent. By a similar calculation, the best-case 
fidelities for FDD and SDD may be expressed as: 



which we could not calculate analyticallySI,. The single- 
cycle results are presented in Table HI an excellent agree- 
ment is seen between analytical and symbolic results, as 
reported in the second and third row, respectively (see 
also Appendix |X] for yet another analytical derivation of 
the coefficient a and the leading decoherence order n for 
CFMG, which agrees well with the results given by the 
other methods). 

An alternative way for determining Fi,(Tc) is provided 
by direct numerical simulation of the total evolution, fol- 
lowed by the process tomography procedure described 
in Sec. IIIDI By restricting the simulations to short r, 
and fitting 1 — Fh{Tc) to a power-law function of r, it 
is then possible to independently determine the leading 
decoherence orders in the ME series. Fig.[l] for instance, 
shows the dependence upon r of F\j{Tc) for N — 5. The 
initial region in Fig. [U where 1 — F^iTc) has a power- 
law dependence on t, indicates the region of validity for 
AHT/ME, the slope of the curve determining the order of 
the best-case decoherence term. The first line of Table U 
shows the leading ME order, rt(ME), obtained in this 
way, which is in excellent agreement with the available 
analytical estimates obtained from AHT/ME. In partic- 
ular, these simulations confirm that Ti!'^^ is the leading 
decoherence term for FCDD2 in the best-case scenario. 



cji_o 



£: -20 




ln(co T) 

^ c 



Fb{T,) 



PDD 



1 - T.^Trs 



PB 



(0) (ah^'^) 



1 - (}T^T^ + 0{T^T^ 



(23) 



FbiT,) 



SDD 



1 - Trs 



Pb(0) [ah 



<4) 



10s(O)) 



(24) 



for suitable real parameters /?, 7. Thus, /c = 2, 4 for FDD 
and SDD, respectively. 

Beside the above analytical determination, the rele- 
vant fidelity decay terms may also be obtained symbol- 
ically - by Taylor-expanding the exact evolution opera- 



tor U{T,) = 



for the appropriate transformed 



Hamiltonians Hj , and collecting coefficients of order r" , 
so that the lowest two contributions to Ti may be iso- 
lated. This approach is required, in particular, to ex- 
tract the best-case leading decoherence order for FCDD2, 



FIG. 1: (Color online) Best-case performance for single-cycle 
DD under the CPMG (plus signs), PDD (circles), SDD (up- 
ward pointing triangles), and PCDD2 (downward pointing tri- 
angles) protocol in the small r region. Bath size N — 5. Solid 
lines are linear fits, with fitting parameters given in Tab. IT] 



TABLE I: Fitting parameters of the best-case cycle perfor- 
mance Fb(Tc) = 1 - K {ujctY'T^, k G R, for a single DD cycle 
under CPMG, PDD, SDD, and PCDD2 in the small r region 
where WcTc ^ 1. Bath size A'^ = 5. Also presented are ME 
predictions by analytical method, n(ME), and symbolic Tay- 
lor expansion, n(Sym). The same agreement holds for N = 3 
and N = 7. 





CPMG 


PDD 


SDD PCDD2 


n(Fit) 


2.00 


4.02 


7.93 11.80 


n(ME) 


2 


4 


8 


n(Sym) 


2 


4 
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FIG. 2: Best-case performance of (a) CPMG, (b) FDD, (c) 
SDD, and (d) PCDD2 starting from an initial state along 
the half-cycle direction {z for (a) and (c); y for (b) and (d), 
respectively). Bath size: iV = 15. r = 0.1, 0.2, 0.3, 0.4, 0.5, 
from top to bottom in all panels (in panel (d), the curves for 
T = 0.1 and 0.2 are indistinguishable). Decoherence freezes 
at sufficiently long time. 



Table U and Fig. [T] clearly demonstrate how DD per- 
formance improves as we go from FDD to SDD and to 
PCDD2. Even though, as already remarked, the CPMG 
protocol [CqZCqZ] does not achieve maximal DD for zero 
bias field, the existence of an approximate integral of mo- 
tion, 5*2, still makes it possible to decouple with high fi- 
delity provided that the initial electron spin state is a S^- 
eigenstate. For two-axis cyclic DD, a preferred direction 
for initialization may still be identified. The FDD pro- 
tocol [CoXCoZCqXCqZ], for instance, conserves Sy in 
the limit r ^ (because Z X ^ Y, which coincides with 
the half-cycle direction); this establishes a Sj^-eigenstate 
as the best-case state for FDD. When several DD cycles 
are implemented, the existence of approximate control- 
induced symmetries is responsible for the decoherence 
freezing phenomenon which we address next. 



2. Long-time behavior 

At long times, none of the above-mentioned analytical 
or symbolic methods is applicable, thus numerical simu- 
lation is required. 

Fig. [2] shows the long-time best-case performance for 
the deterministic control protocols discussed so far. De- 
coherence freezing is clearly seen as a plateau of Fb at 
sufficiently long evolution timesS^!^, the corresponding 
saturation value increasing as the pulse delay r decreases. 
From a control standpoint, decoherence freezing may 
be thought of as signaling the dynamical generation of 
a stable one-dimensional decoherence free subspace via 
j-)]-) 55^56 NMR language, the resulting saturation 

is reminiscent of the "pedestals" seen in the long-time 
magnetization signal under pulsed spin-locking condi- 



-i= -4 



-2-1012 

ln(a t) 

FIG. 3: (Color online) Dependence of the coherence satura- 
tion value on pulse delay for CFMG (plus signs), FDD (cir- 
cles), SDD (upward pointing triangles), and PCDD2 (down- 
ward pointing triangles). Dashed lines are linear fittings at 
small r. 



tive magnetic field is created by the DD pulses, and that 
proper alignment of the initial spin prevents the electron 
from precessing around this direction - thereby experi- 
encing fully frozen nuclear spin fiuctuations. (See also 
Ref. 'Sg for a recent demonstration of a similar-in-spirit 
stabilization effect at zero applied field, leading to a long- 
lived nuclear-spin polaron state via optical pumping with 
polarized light). In practice, decoherence freezing may be 
exploited to optimally preserve a known initial electron 
spin state, through appropriate design of a DD protocol 
with the desired quasi-integral of motion. While a similar 
effect could be achieved by applying a strong static bias 
field, one advantage of DD stabilization is that better 
storage may be ensured by simply rearranging the pulse 
sequence, so that it implements a higher-level protocol 
(from FDD to PCDD2, for example). 

While asymptotic saturation behavior has been re- 
ported for purely dephasing spin-boson models with ar- 
bitrary initial spin statesi^i^i^i^i^, saturation effects 
for a spin bath have only received attention very re- 
cently— i^iSl. Thanks to its inherent simplicity, the 
CPMG protocol makes it possible to gain analytical in- 
sight under the assumption that a simplified coupling 
Hamiltonian is appropriate: 



N 



H 



QSA 
SB 



AS 



k=l 



(25) 



tion! 



,14.57.58 



Physically, one may think that an effec- 



This corresponds to assuming a uniform electron den- 
sity in the QD, allowing one to regard the total nuclear 
spin as a constant. For our purposes, because of the 
large number of spins in the bath, this is also equiva- 
lent to describing the nuclear spin reservoir under the so- 
called Quasi-Static Approximation (QSA), which treats 
the Overhauser nuclear field as a classical random static 
field^. While in principle the QSA is valid for short evo- 
lution time a^^i'^° , it is important to realize that the do- 
main of validity of the QSA is extended in the presence of 
DD pulses, similarly to what happens for an external bias 
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field'^° , and consistent with the fact that the nuclear field 
becomes progressively more static relative to the electron 
dynamics as the electron- nuclear coupling is suppressed. 

By following the steps illustrated in Appendix 1X1 the 
freezing value reported in Ref. IH is obtained: 

^CPMG _ _ 

/ ~ 2r2*2 ' 

for sufficiently small r. Interestingly, the leading power 
of T, n = 2, does not explicitly depend on the number N 
of bath spins. For other protocols in the same range of 
T, numerical results (see Fig. [3|) suggest a similar depen- 
dence of the asymptotic coherence value, 

FfKl-a{aT)", aeM, 

with the relevant values of n being given in Table|TTl Note 
that the characteristic values of r considered here are of 
order of l/cr, that is, uJcT ^ 1 in the simulations, thereby 
well beyond the convergence region of AHT/ME. 

TABLE II: Fitting parameters of the decoherence freezing 
value Ff = 1- a(crr)" for CPMG, FDD, SDD, and PCDDa 
at small r, where o-r < 1 but UcT S> 1. 





CPMG FDD SDD FCDDa 


n(Fit) 


2.0 1.7 2.7 5.3 




ln((n x) 



FIG. 4: (Golor online) Worst-case single-cycle performance 
for FDD (circles), SDD (upward pointing triangles), and 
FCDD2 (downward pointing triangles) in the small r region. 
Bath size: N = 5. Solid lines are linear fittings, with fitting 
parameters given in Tab. IIIII 

case, Fw(Tc) depends on t according to a power-law, 
when T is sufficiently small to ensure that lUcTc < 1 and 
that the AHT/ME approach is vafid. Table Hill f first fine) 
gives the leading decoherence term orders extracted from 
Fig. m which are in excellent agreement with the analyt- 
ical and symbolic predictions (last two lines). Note that 
CPMG, which is not a universal decoupling sequence, is 
not useful in the worst-case scenario, thus we do not con- 
sider it further. 



B. Worst-case performance and arbitrary state 
preservation 

As evidenced by the best-case considerations presented 
above, for sufficiently small r, initial electron spin states 
which are (approximate) eigenstates of the decoupled 
evolution are stable at long times, whereas the spin com- 
ponents perpendicular to the decoherence-free axis are 
lost in the long-time regime. In such a picture, the worst 
case scenario for a given cyclic protocol corresponds to 
initial spin states which are perpendicular to the effective 
half-cycle control axis. Clearly, worst-case performance 
lower-bounds the fidelity of storage achievable for an ar- 
bitrary (possibly unknown) initial state, as required for 
an electron-spin quantum memory. 

1. Short-time behavior 

In order to quantitatively assess worst-case DD fideli- 
ties, we begin, as in the previous Section, by examining 
single-cycle performance. Again, the order of the lead- 
ing decoherence term is determined using three meth- 
ods: (i) exact numerical simulation; (ii) analytical pre- 
dictions based on AHT/ME; and (iii) symbolic Taylor 
expansion. Fig. 2] shows the exact r-dependence of the 
worst-case fidelity for single-cycle DD. Similar to the best 



TABLE III: Fitting parameters of worst-case single-cycle per- 
formance F^{Tc) = 1 - k(ujctYT'^, k G R, for FDD, SDD, 
and FCDD2 in the small r region where uocT < 1. 





FDD 


SDD FCDD2 


n(Fit) 


2.00 


3.99 7.94 


n(ME) 


2 


4 8 


n(Sym) 


2 


4 8 



2. Long-time behavior and randomized DD 

For long evolution times, the worst-case performance 
Fw{T) of six DD protocols (three deterministic and three 
randomized) are summarized in Fig. [5| All schemes lead 
to substantial enhancement of the electron spin coher- 
ence, some of them by more than a factor of 1000, with 
PCDD2 showing the most dramatic improvement^'^. Ac- 
cording to the ME, the leading order term for the FID sig- 
nal is n'-°^ , whereas it is h'^^ for FDD, 7?*^^ for SDD, and 

Ti^^'' for PCDD2, as discussed in Sec. Ill B 21 For random- 
ized protocols, lack of periodicity prevents the definition 
of a time-independent average Hamiltonian, thus AHT 
is not applicable. However, for a given evolution time, 
an effective Hamiltonian and the corresponding leading 
orders to coherence decay may still be defined directly in 
terms of the unitary logical-frame propagato r ^"^i^^ . 
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FIG. 5: (Color online) Worst-case DD performance in the 
logical frame^i^, with r = 0.1. Hamiltonian parameters are: 
tJo = 0, To = and A*' — 15. For deterministic DD, data 
points are acquired at the completion of each cycle, while 
for NRD and FID this is done after every r, for RPD after 
every 4t, and for SRPD after every Sr. Random protocols 
are averaged over 100 control realizations. 



In general, protocols with higher leading order tend 
to give superior performance. However, such a conclu- 
sion does not necessarily hold if a randomized protocol is 
compared to a deterministic one. For example, RPD has 
a lower-order leading term than SDD, so that SDD may 
outperform RPD at short times. However, Fig. [5] shows 
that at long times RPD outperforms SDD, which demon- 
strates the advantage of randomization in suppressing 
coherent error accumulation. The poor performance of 
NRD is expected, since the advantages of this low-level 
DD scheme may emerge only when Q is large, and cyclic 
DD is inefficient. In a closed system^, SRPD has been 
found superior to PCDD2 at long times, but for the QD 
model considered here SRPD does not match PCDD2, 
confirming the fact that irreducible DD groups and slow 
baths are especially favorable for concatenated controlJ^. 
However, it is worth emphasizing that increase in the con- 
catenation level for fixed pulse separation does not nec- 
essarily improve the protocol performance: As shown in 
Ref. 23, PCDD4 may deliver worse fidehty than PCDD2 
if T becomes sufficiently large. A possible explanation 
may be rooted in the DD-induced renormalization of the 
pure-bath terms discussed in Sec. Ill B 21 which for large 
T may become important enough to offset the benefits 
associated with a more elaborated DD cycle. 



IV. REAL-SYSTEM CONSIDERATIONS 

In realistic scenarios, several factors beyond the simpli- 
fied treatment considered thus far will unavoidably affect 
DD performance. Among those related to the underly- 
ing QD model Hamiltonian, the most important are the 
influence of spin-bath dynamics, external bias fields, and 
nuclear spin polarization. We shall address these factors 
one by one, by primarily focusing on the worst-case fi- 
delity of the PCDD2 protocol, which has emerged as the 



best performer for the problem under exam. 



A. Bath size 

In our simulations, the number of bath spins is moder- 
ately large, N < 25, thus it is essential to assess to what 
extent our numerical results might be applicable to real 
QD devices with N ^ lO^'-lO^. In order to verify this, 
PCDD2 simulations have been carried out for baths with 
N varying from 15 to 25, their corresponding spectral 
width being characterized by cr = ^/NA/2 as in Eq. (|20p . 
Fig. [5] illustrates, for different N, the instant of time To. 9 
where Fy^{T) for PCDD2 reaches a threshold value of 0.9 
- as a function of the dimensionless parameter 1/(2ctt). 
It is seen how, by correctly rescaling r (that is, by mea- 
suring T in units of l/cr), the curves corresponding to 
different bath size tend to fall on top of each other, espe- 
cially as TV increases. This provides strong evidence that 
our results should be applicable to realistic mesoscopic 
spin environments upon appropriate parameter scaling. 

It is also interesting to stress that To.g increases ex- 
tremely rapidly as the product ctt decreases and the re- 
gion of very fast DD is entered (note the logarithmic scale 
of the y-axSs in Fig. [6]) . This rapid increase has been an- 
alyzed earlier—. We remark here that a naive extrapola- 
tion from the ME, To g oc 1/r^, which could be expected 
to hold in the limit r ^ 0, strongly disagrees with our 
data in the long-time parameter range of Fig. [HI and that 
a Zeno-type analysis as invoked in Ref. [l^ appears more 
appropriate to explain the observed r-dependence. 



B. Intrabath interaction 

The effect of the internal bath Hamiltonian Hb: 
Eq. ([3]), may become important once the electron co- 
herence time becomes longer than the characteristic 
time scale of the corresponding bath evolution. As re- 
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FIG. 6: (Color online) Bath size effect. To. 9 vs. l/(2(Tr) for 
PCDD2 with different bath sizes A'^ = 15 (black squares), 17 
(blue circles), 20 (green upward pointing triangles), and 25 
(red downward pointing triangles). 
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Time 

FIG. 7: (Color online) Bath dynamics effect for PCDD2 with 
To = (black dashed lines), 0.01 (red solid lines), and 0.1 
(blue crosses). The parameters are = 20 and r — 0.2, 0.25, 
0.3, and 0.4 from top to bottom. 



marked, this effect may be enhanced in principle by the 
PCDD2-induced pure-bath dipolar contribution, given in 
Eq. (fT7|) . In order to assess at which point DD perfor- 
mance start to be significantly affected by nuclear spin 
dynamics, we have performed numerical simulations by 
choosing Tki as uniformly random numbers in [— Fq, Fq], 
and by manually increasing Fq up to values comparable 
to O.lAk, to keep the simulation time reasonably short. 
In this way, in view of Eq. (fT8|) . we cause the effects of 
Hb to be at least a factor of 3 larger than the ones due 

to H for the maximum values of Fq explored. 

The results are summarized in Fig. [71 where a two- 
dimensional 4x5 QD with nearest neighbor-intrabath 
coupling is considered. On one hand, the performance 
of PCDD2 is affected only slightly for sufficiently small r 
and Fq, as expected. Note that for a standard GaAs QD, 
the characteristic time scale for the nuclear spin dynam- 
ics is lOOfis, implying that Fq <C 0.01 in our model. 
Thus, results obtained for Fq = are applicable to typ- 
ical GaAs QDs. On the other hand, as it is also clear 
from Fig. [7l the long-time PCDD2 fidelity deteriorates 
significantly in the presence of a sufficiently fast bath, 
especially for larger pulse delay r. Based on the spin- 
chain results of Ref. 0, randomized protocols such as 
SRPD are capable to be, on average, more robust to the 
effects of the underlying bath evolution. 



C. Magnetic bias field and initial batli polarization 

In practice, a non-zero Zeeman splitting in Eq. ([2|) 
might either be desirable in order to impose a controllable 
quantization axis on the electron spin and/or in any case 
result from residual magnetic fields in the device. As 
numerical results show (see Fig. [Sl[a)), the dependence 
of the PCDD2 worst-case performance is non-monotonic 




FIG. 8: (Color online) Effect of bias magnetic field for differ- 
ent evolution times (top) and effect of initial bath polarization 
for (bottom) on the worst-case performance of PCDD2. Pa- 
rameters are r = 0.3 and A'' = 15. All energies are measured 
in units of A. 

with Bq, by decreasing with small magnetic fields and 
improving for large magnetic fields - this effect being 
more prominent for long evolution times. While this be- 
havior is qualitatively consistent with the emergence of 
a predominantly pure-dephasing process in the presence 
of a strong bias field, the latter also leads, in general, to 
temporal modulation of the fidelity, similar to the effect 
of electron spin echo envelope modulation®"^. 

Interestingly, a non-zero initial polarization of the bath 
spins acts in a similar manner, see Fig. [5Jb). In these 
simulations, we assumed that the initial state of the bath 
is described by a density matrix of the form 

Pb{0) = ^e-^'^^^i^s = [2cosh(/3/2)]^, 
Zp 

where /3 denotes as usual the inverse temperature. Ac- 
cordingly, the bath polarization is defined as the ratio 

1 ^ 

with (•) denoting the expectation over the above bath 
spin state. It has been shown earlier^^ that for the free 
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decoherence dynamics of an electron spin in a QD, a 
non-zero Overhauser field produced by a small nuclear- 
spin polarization is essentially equivalent to an exter- 
nal bias field. Here, a similar equivalence emerges for 
the worst-case performance of PCDD2: the degradation 
trend of the PCDD2 fidelity for small polarizations seen 
in Fig. resembles that occurring for small bias mag- 
netic fields, ajQ < 4. This indicates that the effect of 
the Overhauser field of a weakly polarized bath is also 
equivalent to that of an external bias field in the pres- 
ence of DD pulses. We are unable to further explore this 
equivalence for larger values of the initial polarization, 
due to the fact that exact simulations with a highly po- 
larized spin bath require numerical techniques which are 
beyond our current capabilities-''^. On physical grounds, 
we expect that the degradation trend should stop, and 
improved performance should emerge as p approaches 1 , 
since the fidelity should achieve 100% for a fully polarized 
initial bath state (even in the absence of control pulses, 
for appropriate electron spin alignment). 

V. CONTROL RESOURCES 

We conclude our analysis by addressing the main sim- 
plifying assumptions and requirements implicit in the 
control capabilities invoked so far, with respect to rel- 
evant practical constraints. 



A. Effect of pulse imperfections 

In all simulations presented thus far, control pulses 
have been assumed to exactly achieve a rotation angle 
of (/) = TT in no time - corresponding to zero width and 
infinite strength. In reality, pulses are clearly finite in 
both strength and width, and the implemented rotation 
angle typically deviates from the intended value due to 
both systematic control faults and stochastic parameter 
fluctuations. While a fully accurate error modelling is 
necessarily dependent upon the details of the physical im- 
plementation, our aim in what follows is to gain a sense 
on how stable PCDD2 performance is in the presence of 
different representative control errors. 

The effect of a systematic error in the rotation angle 
may be modelled by assuming that (j) — '"'(1 ~ £)i while 
keeping the pulse instantaneous, and letting e G [0, 1] 
represent the relative flip-angle error. Figure [H] summa- 
rizes numerical results for various error strengths. The 
protocol performance is not very sensitive to this type of 
error, especially for small e (as expected). In fact, the 
electron coherence time is still two orders of magnitude 
longer than Tj* even for e = 0.03, which corresponds to 
~ 5.4°. With phase compensation technique o^^'^^ , DD 
performance might be further improved. 

The effect of finite pulse durations may be accounted 
for by assuming that each pulse is rectangular, with the 
width w and the corresponding strength t:/w adjusted to 




Time 



FIG. 9: Effect of systematic rotation-angle errors on the 
worst-case PCDD2 performance, with r — 0.3 and N = 15. 
The relative flip angle errors are e = 0, 0.001 (indistinguish- 
able from £ = 0), 0.003, 0.01, 0.03, from top to bottom. 

represent an ideal tt pulse for an isolated spin. That is, 
the relevant control Hamiltonian becomes 

H'Jt) = y -(5 • h(t))\h(£T - <) - h(£T - t - w)], 
^ — ' w 

i 

where h{-) is the Heaviside step function. In simulations, 
we have considered the width of the pulse to be up to one 
third of the pulse delay, w < r/3. Figure [TUl summarizes 
results on the worst-case PCDD2 performance for differ- 
ent values of w: Clearly, DD performance depends quite 
sensitively on the pulse width. For a typical QD with 
decoherence time T2* ~ 10 ns, this means that pulse de- 
lays T ~ 4 ns and pulse widths no longer than w = 0.024 
(about 0.3 ns) are required in order to extend T2 by a 
factor of - 100 under PCDD2. 

Improved control design is needed to relax such strin- 
gent limitations, in particular to ensure that high-fidelity 
DD may be achieved with finite bandwidth. While fur- 
ther analysis along these lines is beyond our current 
scopes, preliminary results indicate how an approach 
based on Eulerian controS^^ may allow pulse widths as 
long as T to be employed®^ . A yet more sophisticated ap- 
proach is to resort to numerical pulse shape optimization, 
for instance based on the recently proposed open-system 
gradient ascent pulse engineering algorithm^i^. 

B. Feasibility considerations 

From a practical standpoint, the main requirements for 
DD implementations arc the ability to effect sufficiently 
fast single-spin rotations - along two orthogonal axes in 
an appropriate frame if storage of arbitrary electron spin 
states is sought. While this is a highly non-trivial task, 
single-spin rotations have by now been experimentally 
realized in both gate-defined radiofrequency-controUed 
GaAs double QDs^? and self- assembled singly-charged 
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FIG. 10: Effect of finite pulse widtii on tiie worst-case PCDD2 
performance, witli Tc = 4.8 = 16r + 20w and iV = 15. Tlie 
pulse widths are w = 0, 0.0024, 0.008, 0.024, 0.08, from top 
to bottom. 



(In,Ga)As/GaAs QDs^i^iiZ^ - proposals for further im- 
proving rotation speed and gating time being actively 
investigated in paralleP"'^^. Likewise, multipulse CPMG 
protocols have been successfully implemented on single 
impurity centers in a solid-state matri3*2^i2^i^i2i. In typ- 
ical ESR settings, for instance, magnetic pulses as nar- 
row as 20 ns and gating times of the order of 100 ns are 
currently attainable for typical GaAs QDs at dilution re- 
frigerator temperatures^. 

While the above-mentioned accomplishments and fig- 
ures give hope that full-fledged DD experiments in QDs 
should become accessible in a near future, some addi- 
tional remarks may be useful in connection with the 
prospective relevance of our results to different control 
implementations. Since, as remarked, our main focus 
has been on the zero- field limit {loo = 0), standard ESR 
techniques are not directly viable to effect the intended 
rotations. Rather, the control we are envisioning is based 
on direct magnetic switching, which may be accomplished 
in principle by having access to two independent induc- 
tances oriented along perpendicular axes. In order for 
the required control to be achievable via ESR or coher- 
ent Raman spectroscopy in the presence of a permanent 
static field, pulse widths and separations in the (sub)ns 
range would imply a resonance frequency wq 1 GHz. 
In these conditions, as mentioned, the hyperfine decoher- 
ence process would be largely dephasing-dominated, and 
application of two-axis DD sequences such as PCDD2 
would result in even better performance than at zero 
field - similar to the conclusion reached in Sec. IIV CI 
As also noticed in Ref. 21, a potential advantage of con- 
trol techniques not relying on the presence of an external 
field, however, is applicability over a wider range of pa- 
rameters, along with the possibility to avoiding pertur- 
bations on the spin dynamics in between control opera- 
tions. From this perspective, a more careful feasibility 
study of direct switching schemes in QD devices appears 



well-worth pursuing. 



VI. CONCLUSION 

We have provided an in-depth quantitative analysis of 
the dynamical decoupling problem for a central spin sys- 
tem in the zero- field limit. While our main intended ap- 
plication has been long-time preservation of electron spin 
coherence in a quantum dot by suppression of hyperfine- 
induced coupling, we expect our main methods and con- 
clusions to be relevant for control of nanospin systems 
described by a similar model Hamiltonian. Our main 
conclusions may be summarized as follows. 

For short evolution times and inter-pulse delays r 
which obey the condition uj^t < 1, being the spectral 
cut-off frequency of the total system, analytical results 
based on average Hamiltonian theory and the Magnus 
expansion are in excellent agreement with the results of 
exact numerical simulations. For longer evolution times 
and pulse separations, which are beyond the domain of 
applicability of analytical approaches, DD can still en- 
sure high-fidelity preservation of arbitrary electron spin 
states, as long as typical control time scales are short 
with respect to the spectral width of the total system. 

If knowledge of the initial electron state is available, 
cyclic DD protocols capable to completely freeze electron 
decoherence in the long-time limit may be designed - 
the attainable coherence value depending on both the 
protocol and the pulse delay. 

By studying the effect of important experimental fac- 
tors and control non-idealities - including the effect of 
intrabath interactions, external magnetic fields, and typ- 
ical systematic pulse errors - we conclude that even im- 
perfectly implemented decoupling protocols may still be 
able to significantly prolong the electron coherence time. 
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APPENDIX A: ANALYTICAL STUDY OF 
DECOHERENCE FREEZING UNDER CPMG 
DECOUPLING 

The simplicity of the CPMG {tZtZ) protocol allows 
for an analytical study of the saturation effect discussed 
in Sec. [Tira This wiU be done below via the QSA^SiSO, 
as well as based on a semiclassical random field model. 



1. Quasi-Static Approximation 

Within the QSASii^S, the Hamihonian Eq. ^ reduces 
to Ha Eq. dSll). Let I = Efc Ife' and M = 4. Due to the 
symmetry of Ha and the CPMG protocol, and M + S^ 
are constants of motion. Let |^'(0)) = | t) ® I^j-^) de- 
note the the joint initial state, where the first ket-vector 
corresponds to the state of the electron spin and the sec- 
ond ket-vector denotes the state of the bath. Then the 
evolution of the controlled system only couples the pairs 
of levels I t) ® \I, M) and \ i)®\I,M + 1). (Similar con- 
siderations apply to the case where the electron spin is 
down, initially). 

After n CPMG cycles, the evolution operator is 



U{2nT) 



di d2 

d2 -dl 



where we have defined 



di 
d2 



cos f2- 



t— sm S2 — , 

n 2' 



with C = A^{I - M){I + M +1), B = A{M + 1/2), 
and = . In the best-case scenario, the overall 

system is initialized in the state 

|*(0))b = | T)®|/,M), 

whose survival probability is given by 

\{^m^(2nT))\^ = 1 7 tan^ 61 cos^ 2n6', 



where tan6' = d/Vl — d"^ with d = Im(o?i). 

The initial bath state is a completely mixed state and 
can be rewritten in the basis of |/, M) as 

Pb(0) = ^P(/,Af)|/,M)(/,M|, 

I.M 



with P(/,M) ~ (//DV27rD )e'^ for large N and 
D = Averaging over the nuclear spin states, we 

obtain the survival probability, the input-output fidelity, 
at time T — 2nT, 

r (72 
F(2nT) = 1 - / dIdMP{I, M)— tan^ eicos^ 2ne 

(1 -I- cos4n6') 



C2 



dMAfF(/,Af)^tan" 
B'^ 



For n large enough so that 4n0 ^ 1, we may safely 
neglect the contribution from rapidly oscillating part 
cos AnO. In these conditions, F{2nT) becomes time- 
independent and dccoherence is frozen, 

I r 

F(2nr) ^ Ff = ^ ~ 2 / dIdMPiI,M)^ tan^ 6. 



In the limit of small r, this yields 
1 



f 



1 
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T^A^N = 1 



2^2*2 ■ 



For randomly distributed Ak, A'^N J2k ^1- 

For a single control cycle (n — 1) as considered in 
Sec. nil A H and sufficiently small r so that ^ 1, it is 
straightforward to find that 



FiTc) 



{l~ADT^)e 



-2Dt 



4 1 

5 " 15 
~ 4(1 - i?r2)e-^-V2 

« 1 - + ©(r*^) = 1 - KT^T^ + 0{t*T^), 

consistent with the result reported in the main text. 



2. Classical Random Field Model 

Yet another simple way to describe the saturation as- 
sociated with the CPMG protocol at long times is as 
follows. In situations where the back-action effects from 
the system into the bath may be neglected, a plausible 
approximation is to treat the environment as a classical 
random field. This translates into rewriting the coupling 
Hamiltonian as 

H = B-a , 

where the effects of the randomly fiuctuating field B, 

Bx ~ B cos 6 sin 0, 
By = B smO snY(f), 
Bz — B cos 0, 

are taken into account by averaging over the entire Bloch 
sphere of the system spin, < < vr and Q <9 < 2tt. 

Under CPMG, the unitary evolution after the comple- 
tion of n cycles may be written as 

U{2nT) = (^e'^{-B^X-ByY+B,Z)T^-t(B^X+ByY+B,Z)Ty 



e*^" ) 



where V is the matrix of eigenvectors of U{Tc) and 
exp(±iA), A G M denote the corresponding eigenvalues. 
For a particular initial state aligned with the dominant 
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term of the AHT, that is = B^Z, the fidehty after 
n cycles becomes 



Fb{2nT) 



cos(Bt)2 + sin(Br)2 sin(0)2(l - sm{nX)^) 
cos(Bt)2 + sm{BTf sm{(j) f 



In the long time limit, n = i/(2r) — s- oo, and upon aver- 
aging over and 0, we finally obtain 



Fh^Ff = l 



in agreement with the expression quoted in the main text. 



Notice that the same result may also be obtained by 
restricting the analysis to the two dominant terms in the 

AHT (that is, and by writing the propagator 

after n cycles as 



U{t ^ 2nT) = e 



-i[B^Zt+B^{B^Y-ByX)tT] 



While in principle no a priori reason exists to expect 
such AHT-based description to yield the correct answer, 
a similar treatment successfully describes long-time mag- 
netization "pedestals" in NMR decouplingii. 
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